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ABSTRACT 

BHARAT BHUSHAN PRASAD 
M. Tech. Thesis 

Indian Institute of Technology, Kanpur. 

IMPACT LOADS ON BEAMS ON ELASTIC FOUNDATION 


Modal Analysis and Energy Methods have been used to 
predict the response of beams on elastic foundation subjected to 
impact loads. Experiments have been conducted to verify these 
theoretical predictions/' and they have been found in good agreement. V 
The effect of foundation modulus and external damping on the beam 
response' have been studied in detail. Displacement bounds have 
been computed within and outside elastic range for beams on 
elastic foundation for e certain class of dynamic loading and 
impact loading using Energy Theorems. ■'■- Without going into the \ 

precise details of deformation process, the permanent deformations 

have been obtained for impact loads. The accuracy of these bounds j 

1 

have been checked with a few exact solutions. ^ For a general class {* 
of dynamic loading as well as for impact loads, bounds on displace- 
ment and energy have been computed using variational prinoiples. In 
a few cases, the values of displacement bounds have been compared 




iv 

with the maximum displacement given by Modal Analysis and Energy 
Methods. The wave motion of a beam on elastic foundation in 
bending vibration after impact has been studied to predict how 
the flexural pulse and energy propagate in the beam. It has 
been observed that flexural pulse moves faster in a beam on elastic 
foundation than in a beam without foundation, where as energy 
moves at a slower rate along a beam on elastic foundation than 
along a beam without foundation. 
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CHAPTER I 


INTRODUCTION 

1.1 TRANSVERSE IMPACT ON BEAMS s 

During and after the second World war, scientists and 
engineers have shown an increasing interest in solution of problems 
concerned with the dynamic loading in general and Impact loading 
in particular. The dynamic loading effects on stationary structures 
eg. Railways, Bridges, Runways for landing of Aircraft, Naval 
structures etc., have been the subject of great interest. A lot of 
work has been done on the behaviour of beams, subjected to impact 
loads (Ref. 6 ,7,8, 12, 13, 14, 1 5 » 18, 31, 32, 33), but for beams on elastic 
foundations, very little information is available in . the existing 
literature . 

The duration of impact, the impact force exerted, peak 
values of deflection, stresses, strains, bending moment and shear 
etc. are the quantities of interest for such a situation. In 
addition to other factors, the foundation modulus, the external 
damping etc. play an important role in this soil-structure interaction 
problem. In the immediate neighbourhood of impact point, the 
transient stress distribution is obviously three dimensional? the 
waves are not only propagated along the beam but are also reflected 
from all surfaces. At points farther away, the stresses will be more 
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uniformly distributed, across the width, exhibiting the well known 
bending features with stress variation at a given cross-section 
occuring primarily as a function of depth. This pattern also 
extends to the region of contact point after removal of the impact 
load. 

However, in onler to make the problem amenable to theore- 
tical treatment, it is assumed that the variation of stress across the 
width can be neglected. Ihus finally the problem simplifies to a 
two dimensional problem of plane stress. Additional simplifica- 
tion is brought in by neglecting the rotatory inertia and transverse 
shear. The foundation model is that of Winkler [37] > which 
consists of linear elastic springs. 

1.2 IjITERATUHE REVIEW s 

Lot of work has been reported on impact on beams, 

(Ref. 6,7,8, 12, 13, 14, 15, 18, 31> 32,3?) , but very little informa- 
tion is available in literature for impact on beams on elastic 
foundation. Hoppmann [14,15] used energy method to predict the 
response of a simply supported beam with or without damping 
resting on Winkler model of foundation. He observed that with 
increase in the modulus of the elastic foundation supporting the 
beam reduces the amount of energy absorbed by the structure. He 
observed that this is due to the fact that the fundamental period 
of vibration of the structure becomes small in comparision with 
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time of impact. Impacts by relatively small masses moving with 
low velocities develope bending strain of practical concern. 

For cases of impact on simply supported beans on elastic 
foundation, in which several modes are excited, the maximum defle- 
ction rapidly decreases with increase of foundation modulus whereas 
the strains decrease at a much smaller rate. If the fundamental 
mode of vibration of the beam is the only one appreciably excited 
by the impact the strain is proportional to the deflection and 
both decrease at the same rate vdth increasing stiffness of founda- 
tion [ 14 1 • 

External damping has a rapidly decreasing effect on reducing 
deflection and strain, as the number of the mode increases [14] . 

A very large portion of the deflection is due to the fundamental 
mode whereas a considerable portion of the strain is due to the 
higher modes. The co -efficient of restitution for the mass striking 
the beam increases very slowly with increase of foundation stiffness, 
if several modes of vibration are present. 

1.J SCOPE OP THE PRESENT INVEST! CATION s 

Using Model Analysis and Energy Methods, the responses of 
Euler -Bernoulli beam resting on .Winkler model of foundation to impact 
loads have been analysed. Chapter 2 deals vdth analytical investi- 
gations. The effect of external damping, foundation modulus on beam 
response have been studied. Experiments have been conducted to 
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verify the theoretical results. Chapter 5 deals with experiment al 
setup and experimental measurements. The Euler-Bernou 111 -beam on 
Winkler model has been similulated in an experimental model consis- 
ting of an uniform mild steel beam and fairly linear uniform elastic 
springs . 

For a certain class of dynamic loading including impact 
type of loading , displacement bounds for beams on elastic founda- 
tion have been computed using Energy Principles. The analysis of 
impact problems finally boils down to evaluation of initial velocity 
distribution or initial Kinetic Energy input. For a given initial 
velocity distribution, the displacement bounds within and outside 
elastic range have been compited in Chapter 3* These values have 
been compared with exact solutions for checking the accuracy of 
these bounds. The impact load is also characterized as a dynamic 
load q(x, t) using Hertz law of impact. The variational principles 
have been used in Chapter 3 to compute the ui splaceraent bounds of 
an Euler -Bernoulli beam resting on Winkler model of foundation, 
subjected to a general class of dynamic loading q.(x,t) including 
impact type of loading. 

How the flexural waves x^op agate' after the impact and how 
the energy propagates along a uniform beam on elastic foundation 
have been studied in Chapter 4 . It has been shorn that the energy 
propagates with a smaller velocity along beams on elastic foundation, 
as compared to beams without foundation. The general conclusions have 
been drawn in Chapter 6 along with scope for further research. 



CHAPTER 2 


ANALYTICAL INYESH GAIT ON S 

2.1. MODAL ANALYSIS (FREE FREE BEAM ON ELASTIC FOUNDATION WITHOUT 
D PIPING ). 

In this analysis Beams on elastic foundation are considered 
without external damping. A method is developed wherein the 
vibration equation of beam on elastic foundation is solved and the 
applied impulse is replaced by a suitable initial condition. The free 
free beam resting on an elastic foundation is treated as undamped one 
dimensional continuous system. The vibration of such a system is 
represented by the partial differential equation 

a 2 w(P,t) 

£ [W(p,t) ]+ MCP) — = o (2.1.1) 

at 1 ' 

where Lis a linear homogenous self-adjoint differential operator 
consisting of derivatives with respect to the spatial co-ordinate P, 
but not with respect to time, and l(p) is the mass distribution 
function. The above equation reduces to 


3x 


3 W(x,t) 2 V( t) 

[ El _ - ] + K(x) W(x,t) + p A(x) - — = 0 


Cx) 9x 


2 


at/ 


where the linear differential operator £ has the form 



S 2 

— ] + K(x) 
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and M(P) = PA(x). 

For a uniform beam, on a uniform elastic foundation the above 
equation reduces to , 


„ T 3 4 W(x,t) 2 

EI 1 4 — + PA + KW(x, t) = 0 


( 2 . 1 . 2 ) 


Here xt can be shorn that the operator £ is only positive and not 
positive definite, whereas Mis positive definite making it a 

semi-definite system. The solution of the above equation can be 
written as 


)=I X rx).n (t) 


ijaci separating the variables in (2. 1.2) 


A 

d \ 



4 

a X = 0 

n n 


and 


T/h ere 




2 /a 2 ) - (K/RI) 


a 2 = El/p .A 

X n = Eigen function 

w n ■= Eigen frequency 
= Eigen value . 


(2.1.3) 

( 2 . 1 . 4 ) 

(2.1.5) 

( 2 . 1 .6 ) 



Solution of Equation ( 2 . 14 ) yields, 

= C . e , where s = + a n and s = + i 
satisfying the boundary conditions for free-free ends as 


? 


FI 


3 W 




x=o 


x=o 



x=L 


El 



= 0 


x=L 


(2.1.7) 


Then the frequency equation becomes 


Cos a L. Cosh a L = 1 
n n 

This transcendental equation has been solved numerically for the 
eigen values 


0 0 

1 0 

2 4 . 73 OC 4 QB 

3 7.8532C 46 

4 IO. 996 O 78 

5 14.1371655 

6 17 . 27 8 596 


ind for higher values of n, approximately for n > 6 
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and w n = ^ fl ’V+ CK/EI) ( 2 . 1 . 10 ) 

Thus it can "be observed that frequency equation ( 2 . 1.8) yields a 
zero double root. 


ci = 0 and a = 0 
° 1 


for a n = 0 , equation ( 2 . 1 4 ) reduces to 


<i 4 \(x) 


m 0 


dx 


which has the general solution 


X 


n 


D^+ Dg'x + D,x^ + D. x^ 


( 2 . 1 . 11 ) 


( 2 . 1 . 12 ) 


It is easy to see that to satisfy the boundary conditions given by 
equation (2.1. 7 ), I> 7 = D, = 0. Hence corresponding to zero double 

J nr 

roots, there are two modes 

X Q (x) = A q (2.1.15) 

X 1 (x) = iL, ( 39 / L - ■§•) (2.1.14) 

where X Q (x) is selected as symmetric and X^(x) as an anti -symmetric 
mode, as this particular linear combination renders the eigen function 
X Q (x) and X^(x) orthogonal. They correspondingly represent rigid body 
transverse translation and rigid body rotation about the centre of 
mass of the beam. Since ® =a Q = 0 is a double root, any combination 
of X q (x) and X^(x) is also a natural mode. This explains why it is 
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possible to choose X^ (x) to represent rotation about the centre of 
mass rather than any other point. Such unrestrained beams are 
usually treated as semi-definite system [23 ] . 

The frequency equation (2.1.8) further yields an infinite 
sequence of discrete eigenvalues, a ( n = 2,5, ...<») which corres- 
pond to the eigen-function. 


X*(x) = [ (Cos a n L- Cosh a n L) (Sin a„x + Sinh a^x) 


(2.1.15) 

- (Sin a n L- Sinh a^L) (Cos a p x + Cosh a^x) ] 

It is interesting to note that eigen-values and eigen -functions are 
the same as for a free-free beam without elastic foundation. On 


the other hand the natural frequencies are given by equation (2.1.10) 
which are different from natural frequencies of the same beam without 
the elastic foundation. The effect of modulus of foundation dimini- 
shes as the mode number increases. Thus X (x) = X + X„+ X *(x) . 

n v ' o 1 n v/ 

Here X^(x) satisfies the boundary condition as well as the governing 
differential equation. As this semi -definite system has the mass 
distribution function 'M(P) = Pi(x) = m(x) = m which is positive 
definite, and as this is a self-adjoint problem, and X (x) is also an 
admissible function then [23J 


/ m(p) ^(P) X k (P) dD(P) = £ nk 
D 

where 6 is a Kronecker delta 

D denotes the domain. 
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for n , k, 8nt - 1 

n / k, 4nk , 0 


and Xy. is orthogonal to X 


n 


) 


m 


X n (x) ^(x) dx = <5 


nk 


For X (x) = X = A 
n v o o 


m A . A . dx 
o o 


= 1 


A = + 


4 l 


Again for X (x) = X, = A,, (~ - {La 
n 1 1 v l 2 


then 


^ Hi A^( - -|-) 2 dx = 1 


n.A 2 .(L/l2) = 1 


*1 


/_12 

mL 


( 2 . 1 . 16 ) 


( 2 . 1 . 17 ) 


Similarly fox elastic inodes, X* = F n . X^fx) 


L 

/ m(X n *) 2 dx =• 1 
0 



1 _____ 

(Sin a L • Sinh a 1) 
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X n ( x ) = Z=T +/ ~Sr( f -2) + — f (° os a n x + cosh “ x)- 

v mL /mL 

'Cos a L- coda a L 

/ ■ II XI \/ \ 

V A sin. cvx+Sinha x) 

sin a i-sinh a n 

n n 

. ... (2. 1.18) 

Eq.(2.1.5) yields the following solution, 


«n(t) ■ °n 003 V + \ Sin r, n * 

CO 

So w(ic,t) = I X (x) {C cos w t + D Sin w t } (2.I.19) 

n=2 n n n n n 

With initial condition 
W(x,0) = 0 


H = 0 for t o 0, x f L/2 (2.1.20) 

If - V o for t . 0, x - 1/2 

Prom Eq.. (2.1.20) the constant = 0 
So Equation (2.1.19) reduces to 5 

CO 

W(x,t) = Jo D n sin V (2.1.21) 

Defining G = w . D 
n n n 

00 

W( x, t ) « l — - . G X (x) Sin w t 

n=2 n n 


( 2 . 1 . 22 ) 
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The evaluation of the constant G actually involves the replacement 
of the impulsive forces hy an appropriate initial velocity condition 
for the beam [8 ] . If the striker is considered to impart a velocity 
indi stingui shah le from v q . to an infini tesimal section of the beam 
just under the contact point, then the momentum of tills section and 
the striker should equal the initial momentum u t , that is 

/ <x, 0> dQ = m 2 ^ o ( 2 . 1 . 25 ) 


where integration dQ is carried out with respect to the total mass 
of the striker and beam 


aw 

5 t 


(*,0) - l & n X n = *(ac) 


(2.1.24) 


n=2 


Multiplying both sides by another eigen function X^. 


1 G n-\- x k- *«• \ 

n=2 


(2.1.25) 


If Eq,.(2.1.25) is multiplied by m 2 and evaluated at X = L/2 and is 
added to the integral of Eq,.(2.,1 .17 ) with respect to the beam mass 


dm^ = (m^/L)dx, then in view of symmetry, the resulting expression 
can be obtained as 


L/2 

CO “ 

t \ £ 2 / .X n .X k (m 1 /i) dx + m 2 X (1/2) \<&2) ] 

n=2 

° L 

2m 

= -— / 4(x) 1^. dx + m 2 .tf;(L/2) X k (l/2) (2.1.26) 


0 



Eg.. (2.I.I9) represents really an integration with respect to the 

\ 

total mass as suggested "by the Iijulse-moiB ntum equation (2.1.15). 

For the orthogonal modes, it can be shown that left hand side of 

Equation (2. 1.26) vanishes for n f k, and hence 

L 



m l 

— J i[»(x) X n dx + m 2 X n (L/2) ^(L/2) 

2 (2.1.27) 


L 

p 

-I \ dx+ “2 VV2) \(V2) 

0 


Using the initial condition given by Equation (2.I.14), 


G = 
n 


m 2 V o 






(L/2) 


L 

/ 


(2.1.28) 


X dx + in X 


n 


xi 


( l/s) 


and 



Cos’" a 

L. 7 — - (Cosh (v I- Cos a l), ra » pil 

Sin 4 a L ^ n 1 

n 


“ 1 X n (x). X n (L/2).Sin w^t 

n=2 n Cos 2 a L(Cosh a L-Codi a L) +m. 0 7? (1/2) 

pAL r— = n n Z 11 


>/(>:, t) = m 2 v 0 
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l tf(x,t) = mv 


l ± 


X (x). X (1/ 2) Sin w t 


'iO . w _ 2 _ 

n =2 n Co s a I> 

pAL n 


— — (Cosh a 1-Cos a l)+ m X (L/ 2 ) 

Sin 4 a 1 n n 2 n 

n 

( 2 . 1 . 29 ) 


The impact force F(t) is evaluated "by knowing the accelera- 
tion from equation ( 2 . 1 . 2 1 ) , 


3 W / / \ 

m 2 • — 2 * ( V2 » 

3t 


(2.1.50) 


Bending Moment is given by 


I( x, t ) 


- EE(o n r W(x,t) 


(2.1.51) 


Sheer force is given by 


S(x,t) = - El 9 


5 W(x,t) 


J~ “ " El (<£) W(x,t) 


( 2 .I. 52 ) 


Bending strain is given by 


e (x,t) = 


H 

2 2 


(2.1.53) 


where H = depth of beam. 


Hie stress is given by 

a (x,t) = E. H/2. 


(2.1.54) 
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2.2 MODAL ANALYSIS (FEES- FEEL BEAM WITH DAMPING) 

Hie equation of motion for the transverse vibration of 
beam on elastic foundation with external damping, can be expressed 
as 


[ El 


9x 


(x) 


a 2 w(x,t; 


2 

] + K(x) W(x,t) + P(x) A 


9x 


rt 

9 W(x, t) 

0 —— “ 0 


where G = Damping factor. 

For a uniform beam on a uniform homogeneous elastic foundation, 
Eq. (2.2.l) reduces to 


( 2 . 2 . 1 ) 


El 


IjXsO) 

9x r 


+ C 


9W(x,t) 


3t 


+ PA 


:i(.x>t) 

9t 2 


+ KW( x, t) = 0 (2.2.2) 


Separating the variables in Eq.(2.2.2) 


4 

d’X 


* 


dx 


n\ 


0 


( 2 . 2 . 3 ) 


and 


1 > 

dt 2 


C dq n ^ 2 „ 

+ PA dt + 'n “ 


(2.2.4) 


The solution of ( 2 . 2 . 3 ) will be same as Equations (2.1. 7 ) and (2.1.8). 


Putting e = -p^ , Equation ( 2 . 2 . 9 ) will yield two roots given by 
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„ l , 7 2 2 

CL =^{-0+^c -• 4 W -l 

■ **** Y"i J 

1,2 


Further introducing 

c 

£ = r — = damping co -efficient 
2 w n 

then 1 1>2 > - t\ ± A 2 -1 w n 

Hence three conditions can arise 


(2.2.5) 


(2.2.6) 


(a) £ > 1 , making the system over damped 

(b) £ = 1 , the critical!:’ damped System. 

(c) 5 < 1 , undsrdamped system. 


2.2(a) O verdamped System s 

As 5 exceeds unity, the two taa+s remain on the real 
axis and are unequal. The equation for i becomes 


vV W^-i) VtV .(-e - s? -i)k„ * 


M < 4 * t O / r / 

or Si = ® 5 n ( c n Gosh -1 w a t + ! n Sinh w t) (2.2.7a) 

Equation ( 2 .2 ,7a) represents aperiodic motion. So solution of ( 2 . 24 ) 
can bo written as 

°° 

W(x,t) l X (x) e n (C Cosh /£ 2 -1 w t + I) Sinb/^-l w t) 

XI XX lx XX XI. 

xIkc v 

.. .(2.2.8a) 
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Using the initial condition (2.1,20), 'evaluation of C 

n 

and D p is done exactly in the same way as done for the case without 
damping . Finally, the solutions can be written as 




Shear force is given by 


S(x,t) = - EE 




Bending strain is given by 



( 2 . 2 . 10 ) 


(2.2.11 1 


(2.2.12a) 


(2.2.13a) 
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strewn, o (x,t) . E . H/2 . 9 2 W(x,t) 

2 


( 2 . 2 . 14 a) 


2.2 (To) UHDE ED 1 MEED SYSTEM ; 

As 5 is less than unity, the two roots are complex 
conjugates, the solution for becomes, 


Si “ ( \ 


+iw 
e n 



+ B 


n 


~i 

e 


w 

n 


/i~^t) 


e 


- F-w t 

T1 


or 


a 

n 


V (G n Cos V + \ Sin i t) 


where v, 

d 


w n 


✓ 1 



(2.2.7b) 


So solution of Equation (2.2.4) can be written as 
\ (( - 5 w t 

AX,t ~ n=2 Xn ^^ S ^ C n Cos w d t + D n Sin w d t) (2.2.8b) 

Using the initial condition given by Equations (2.1. 13) and ( 2 . 1 . 14 ), 

evaluation ofC n and ^ is done exactly in the same way as done for 
overleaped system. So 


¥(x,t) =: 


v 

2 0 


I 

n=2 


„ . “?w t 

X n (x).X n (1/2) e n 


w. 


Sin w,t- 
d 


PAL Cos « n I -(Cosh oL- Cos a L) 


Sin^'a^L 


+ m 2 X n( L /2) 


(2.2.9b) 
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Hie impact force is evaluated "by knowing acceleration 


from equation (2.2.9b) 


P(t) =- m Ul (L/2, t). . 

^ o 
9t 


(2.2.10b) 


Bending moment i s given By 


M(x, t) * - EE 


9 V x. 


(2.2.11b) 


Shear force is given by 


S(x,t) = - EE 


3^W(x,t 


(2.2.12b) 


Bending strain is given by 


(x,t) =. 


H 9 2 W( x. 


(2.2.13b) 


Ut) .s.ijfei 


(2.2.14b) 


2 . 2 . ( c ) CRITIC ALLY D A&/IEED SYSTEM 

If the damping coefficient 5 equals unity, it corresponds 
to critical damping. The solution for q becomes 


o = (A + B t) 
hi v n n ' 


(2.2.7c) 


Hie final solution can be written as 


w(x,t) = m c v I ~ 


(L/2) .t. e 


2 o w 2 / . p 

n=2 n PAL Cos (Xn L(Coah c^I* - Cos a/L) + (L/2) 


Sin 4 a L 
n 


(2.2.9c) 
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Hie impact fome is evaluated by knowing the acceleration 
from equation (2.2.9c) 


F(t) = - m rw (1/2, t) 

^ 2 
9t 


(2.2. I0e) 


Bending Moment is given by 

2 


M(x,t) = - El 


Sheer force is given by 


Ox 


(2.2.11c) 


S(x, t) = -El 




Ox' 


.3 


(2.2.12c) 


Bending strain is given by 


e(x,t) = 


stress, o (x,t) = s. H/2. 


H 3 W(x,t) 


Ore 


2* 

VV 




X. 


(2.2,1 3 c.) 


(2.2.14a) 


2.3 MODAL ANALYSIS (SIMPLY SUPPORTED BRASS) 

(without damping) 

Ihe equation of motion for the transverse vibration is 
same as given in Equations (2.1.1) and (2.1. 2 ) with the boundary 
condition 

W = 0 , x = 0 

¥ = 0 , x = L 


(2.3.1) 


-El _0_W 
9x 2 


- El 


36=0 




Ox 


( 2 • ‘j . 2 ) 


Xr= L 



21 


As obtained in the Section (2-2), the solutions are 




Sin 


nrrx 

L 



and 


. - 

n l pA j 


• n\ 4 K 1/2 
' L 4 + EI 


(2.3.3) 

(2.34) 


(2.3.5) 


^ - °n Gos w n t + D n Sin w n ( 2 . 3 . 6 ) 

So 

¥(x,t) = Sm (C ft Gos * n t + Sin w n t) ( 2 . 3 . 7 ) 

¥ith thS help 0f initial condition, C n ,D n are evaluated exactly in 
the same may as for Free-Free case. The final solution is 


W(x,t) = 


m 2 v 


o l ^ Sin 21* 


Sin wr/2 


n=l w n 


pAL 2 

— + m 2 Sln n7r/2 


(2.3.8) 


MODAL ANALYSIS (SIMPLY SUPPORTER WITH DAMPING) 

( 2 . 4 a) For € > 1 making the system oveidamped, referring to Eq,.(2.2.4) 
■the solution for c[ "becomes 


-?w n t 


q * 1 = 6 CC n CoSh /52-1 V + D n Sinh ^ 2 “1 v^t) ( 2 . 4 . 1 a) 
With evaluation of constants, as for Fzee-Free case 


uu 

W(x,t) = m 2 v Q l 


Sin £££ Sin 
L 


w 

n=l n 


__ “?W_t 

xiTr xi rrr — *“* 

T“ e • Sinh /^“l- w t 
_£ n 


pAL A m c- 2 

*X" m 2 Sln nTr /2 


(2.4.2a) 
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(2.4b) For 5 < 1 , making the system underdamped, the solution 
for becomes 

~ ?' T ji ^ 

q n = e ( c n Cos w d t + Djj Sin w d t) (2.4.1b) 


where w, = w /1-r 2 
an * 


With evaluation of constants, as for free-free case 


00 , 

W(x,t) = m 2 v 0 l i 


~Ew t 

nnx c . nir ^ n . 
bm — — Sm — e Sm w, t 

L 2 d 


n=l d 


AL . 2 , 

— + Jn 2 Sm n-ir/2 


(2.4.2b) 


The impact force F(t) is evaluated by knowing the 
acceleration from Eq. (2.4.4) 

Bending Moment is given by 


-El 


3 _W (x,t) 
3x 2 


-EI(o n )- W(x,t) 


(2,4.3b) 


Shear Force is given by 


-El 


r%(x,t) 


9x 


3 = - EI(« n )° W(x,t) 


(2.4.4b) 
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(2.4c) For 5=1, making the system critically damped, the solution 
for becomes 

*“W 

Si = CAn + B n t) e n ( 24 . 1 °) 

With evaluation of constants, as for Free-Free case 


W(x,t) = m2 v 0 l 


c . HTTX c . H7T 

Sm — ~ Sm — r 
__ L 2 


t* e 


-w t 

n 


n=l v n 


pAL/2 + m 2 Sin nrr/2 
2.5 ENERGY METHOD (FREE FREE BEAM) 


(2.4. 2o) 


The impact is characterized here as time dependent load 
designated as F(x,t) and the governing differential equation for 
a uniform Free-Free beam with external damping takes the form 


El 


a 4 wfx. tl 

. 2 
3x 


+ p A 


rw(x,t) 

9t 2 


+ C 


3W(x,t) 

at 


+ KW(x,t)*F(x > t) (2.5.1) 


And if impact happens to be in the middle of the beam, F(x,t) = 

F^ j t), then the term F(x,t) may be replaced by an equivalent point 
force which depends only on time. 

As done earlier, the cross-sectional dimension . of the beam 
is considered small in comparison with the length and the effect of 
shear deformation and rotatory inertia is neglected. The boundary 
conditions are 




= 0 


x=0,L 


x=0,L 


( 2 . 5 . 2 ) 
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The initial conditions are expressed as 


W( x, 0) . 0 

awl 


(2.5.3) 


Neglecting contact deformation, vibration of the striker and the 
static deflection due to attachment of mass nr,* the solution of 
Eq. (2.5.l) is given by 


W(x,t) = l X n (x) q n (t) (2.54) 

n=° 

where X n (x) are eigenfunctions satisfying Eq. (2,5'.l) and Eq . (2.5*2), 
as obtained earlier in Eq. (2. 1.18). 

Here q^( t) is determined by substituting in the Lagrange's 
equation the value of Kinetic Energy T and potential energy Y of 
the beam-foundation system. The Lagrange's equation is written as 


i- E 3I^ _ 

dt 1 * J 


9T 


9^ 9q n Sq n 

where D = Energy dissipated due to damping. (Dissipation function) 


3V 3D 


1 


(2.5.5) 


V-ir ma 


— (y ) 
d n ; 




Expressing 

1 00 , L ~ 

T = J l tq n ) PA / (X n ) Z dx 

n=o 0 

V = V BEM + V P0UI©ATI0N 


(2.5.6) 
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+ I K f (q n ) 2 / (X n ) 2 dx 


n=0 


(2.5.7) 


d - 1 L c ^ n ) 2 Av 2 d * 


n=0 


(2.5.8) 


in the Lagrange’s equation (2.5*5)» the resulting equation is 

/ L (X") 2 dx 

q + * f C ] rEI _0_ K 

H n 4 n l pA J L c 


PA ; L 2 

/ v dx 
0 




F(t) X n CL/2) 

pA / L X 2 dx 
0 n 


(2.5.9) 


Using Laplace Transform and with initial condition given by Eq, (2.5.5) 


W(x,t) = 


1 ~ f c -qct-T) 

* Jo 0 F<I)e 


(2.5.10) 


0 


Sin 6 n (t-T)dT * 


where 


L Cos 2 cx_ L 

/ X 2 dx = L 


0 


Sin 4 a L 
n 


(Cosh a n L - Cos c^L) 


6 .!£- 
n 2 pA 

e = [w 2 - i 2 ] 1 '' 2 
n n n 


(2.5.11) 


(2.5.12) 
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j ^ e Sin 0 n (t-r) dr is Luhamel' s (Convolution.) 

0 


integral whose evaluation depends upon F(t). So it is necessary 
to develop an analytical expression for the contact force F(t) in 
in terms of the known physical parameters, so that this Duhamel* s 
integral may he evaluated. Based on Hertz [10] , Law of impact, 
impact forces have been given by Lee [18 ] , Zeener and 

Feshback [59] . Selecting 

F(t) =■ n> 2 v Q (l+e) — Sin 0 < t < 

L L 


t > T 


(2.5.15) 


where = period of contact and e = co-efficient of restitution. 

W(x,t) can be evaluated by putting F(t) in Eq.. (2.5*10) and performing 
the integration. But still the deflection W( x, t) depends upon the 
unknown co-efficient of restitution ' e' . Formula should 
be known so that 1 e* may be determined in terms of known physical 
parameters. This has been evaluated in the same way as suggested 
by Hoppmann [14] for simply supported beam. The contact period 
has been computed as suggested by Hertz and Raman [ 28 ] 


2 c i 2 > 2 

Cl-V ) ,/<- 

3.28 ] 1/5 


V r 


(2.5.14) 


where r = radius of the striker. Finally, 


Trm 2 v (1+e) « \ CL/2) Sin 4 a n L.I n .X n 


E ll u 11 I* , , 

— (2.5.15} 

n=0 Cos z a^L (Cosh a n L-Cos o^L) 


W(x.t) = 


2p AL T l 
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For t < T_ 


-v 


1 = 

EL 


? 2 2 47f2 6 n 

r w 2 _ %-) 2 + 

L n T l J t 2 


2 5 TT ^ 0 7 TT^ 

U-f- Cos V . — (2<-< - ^)Sin e n t] 
L Tj n I t 


2 ^ it . „ 2 

[ Cos I* + (w 2 - -j) Sin It] } 

T L T L n T *L 

L 


(2.5.16) 


For t > T t 


“V 


X n = 


, 2 tt 2 4it 6 n 

(w - — — + — 

n _2 -,2 

1 L *L 


2 ,2 COS e n t <■ J n ) 


6 T l S tt 

p n =[ 2 e n - 2-Cos0 n T r 

n T t L T t 0 


n T L - - 

TT 0 0790 ^ ^ 

'n* l T— C2i n - »; + "‘ /T F Sin 6 n T L * — ) 


L n 


J = 
n 


V r L 

TT0 2 °tJ^L ^ ^ TT 

t~ < 25 n-”n* 3>- 005 6 nV 2a ’ 9 nV TTe - 

t. e t: l Ln 

L n L 


R ] Sin 6 T t 
n J n L 


\ ■ ei - “n + '7 t 2 ) 


(2.5.17J 
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The Bending Moment is given by 


M(x,t) = -El 


3 2 W(x,t) 


3x 


Shear Force is given by 


S(x,t) = -El 


gfwcx^t) 


3x~ 


Bending Strain is given by 


, - H 3 2 W(x,t) 

eCx,t) = — ^ 


3x 


The stress is given by 


0(x , t) . e « . ffeil 


3x 


(2.5.1s) 


(2.5.19) 


( 2 . 5 . 20 '' 


(2.5.21) 
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2.6 RESULTS m DISCUSSION. 

Numerical Results have been obtained for the cases discussed 
in sections (2.1 to 2.5) and are presented in Figs. (2.1 to 2.1 0) which 
show the variation of Deflection, Strain, Impact force, Stress 
with velocity of impact, v q and with time. The summation of 
series have been done for first ten harmonics using IBM The 

convergence of series, Eg.. (2.1 .29), Eg.. (2.2.9a), Eg.. (2.2.9b), 

Eg.. (2.2.9c) for free-free beam on elastic foundation was very good. 
However the convergence of series expressed by Eg.. (2.3.8), Eg.. (2.4.2a) 
Eg.. (2.4.2b), Eg.. (2.4.2c) for simply supported Beam on Elastic 
Foundation was not that good. 

The effect of foundation modulus is to increase the natural 
freguencies of vibration of the beam, however this effect diminishes 
as the mode number increases. The natural freguencies of the free-free 
beam on elastic foundation has been listed in Table (5.1). 

Increment in impact velocity, v q , increases the various 
quantities, like, impact force, deflection, strain, stress etc. as 
shown in Fig. 2.1 to 2.10. 

The effect of external damping has been shown in Fig. 2.1 to 
2.10. In case of £ >1, the motion is no longer oscillatory, being 
exponentially decreasing with time, making the motion aperiodic. 

In case of g < 1, the motion is oscillatory with diminishing 
amplitude. The beam response is greatly reduced by the use of 
external damping. 
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The beam responses, (Deflection, Strain, Stress etc.) as predicted 
by Model inalysis are higher than that obtained by Energy method for 
Free-Free beam as well as simply supported beam as shown in Fig. (2.1 to 
2.10). The beam response for simply supported beam by Modal Analysis 
has been compared with Hoppmannfs results (Energy Method) in Fig. (2.9). 

Here also the results by Modal .Analysis are higher than by the Energy 
Method. 
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CHAPTER 3 


DISPLACEMENT BOUNDS POE BEAMS ON ELASTIC FOUNDATIONS 
3.1 INTRODUCTION 

Prom the viewpoint of Engineering Analysis, and design, a 
principal feature of the response to transient loading or impact loading 
is the permanent deformation that occurs in the Beam on elastic 
foundation. It is therefore of general engineering interest to have 
techniques at hand which establish hounds on the overall deformation 
in dynamic plasticity problems, but which, at the same tine allow 
some or all of the complicating details of the complete response to 
be ignored [29] . 

Even for deformations within the elastic range, the bounds 
on displacement and energy, is of great help in doing the optimal 
design of such beams on elastic foundation, whereas it may not be easy 
to keep these bounds as the objective function and do the conventional 
optimal design^ Namely a special design problem of this kind, may be 
stated as follows : for a given bound on deflection under specified 
dynamic loading, what design provides minimum total weight ? [1, 16, 38] 
The design which minimizes this response bound for dynamic load 
has been shown to be quite different from the optimal design for 
minimum deflection under static loading [26] . Here an attempt 
has been made to compute upperbound on deflection of a beam on elastic 
foundation by means of an energy approach and some suitable inequalitie 
The energy has been taken as functional, and for obtaining a bound on 
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displacement of the structure, an inequality relating the functional 
to the displacement is required. Since the actual motion is not known, 
this inequality must hold for all possible displacements. The 
determination of the optinum inequality of this type is given in 
section (5.%) As the objective of this chapter is to find the upper- 
bounds on deflection of the beam-foundation system at the interface, 
damping has been neglected (as damping is known to reduce the amplitudes) 
in the present analysis. 
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3.2 MARTIN* S APPROACH : 

(■ 

For Hearns resting on an elastic foundation subjected to 
dynamic loading it is generally not possible to -write a simple 
analytical expression for the maximum displacement at a given point 
on the beam. But for a certain class of dynamic loading, it is 
possible, however to formulate a comparatively simple analytical 
expression for an upper bound on the maximum displacement at a 
point. For the elastic structures, such imethode have been ex:amined 
by Martin ['19, 20,21,22] . Martin has illustrated the displace- 

ment bound principles for elastic and inelastic continuua, which 
are in fact inequalities relating equilibrium stress fields and 
compatible displacement or velocity field in the continuum. . 
Treating this beam as a continuum, considering a stress field . . 

jL^J 

such that 

3c 8 

+^=0 (3.2.1) 


at each point in the body and 
0 ij a j = T i 


( 3 . 2 . 2 ) 


at each point on the surface. The stresses c are thus in 

internal equilibrium with body forces F^ and in external 

s 

equilibrium with surface Traction T^ , where 1 a^ ’ is a unit 

outward vector normal to -the surface at the point under considera- 

s 

tion. Further a. . at ea*h point should not ‘violate the initial 

1 J 
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yield condition. Then considering a continuous displacement field 
u i 311 d strain field , such that at each point 





3Xj_ 


(3.2.3) 


Further it is required that the displacement and strain, history 


following a path from e = 0 , u = 0 to s. , u. ° should t 

1 ii i ini 


13 1 


icio tin . 


, S S C G 

Since a . . , T. , F. .ore in equilibrium and e. . , u. 


13 ' 1 1 


13 1 


s,re compatible, by virtual work principle 


/T 5 u. C dA + / F. u, C d¥ = / 0 ., 5 e . 

A 1 1 xr 1 1 -V 3 13 


(3.2.4) 


Por the dynamic problems under consideration, the conditions ares 


1 . At time t < t Q body forces F^ act on the continuum, 

causing displacements and strains, u. 1 , e .. ! such 

1 13 

that all stresses in tlii s state are within yield 
condition. 


2. At time t' = t Q , the continuum accquires velocity u. 0 . 

3. For time t > t , the surface tractions are zero over 
part of the surface and the displacement rates are 
zero over the remainder of the surface. 
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4° for t > t the displacement, velocities and strains 

"fc * "fc "fc 
are given "by m , m, 

Then the energy balance equation relating the states at 


t and t, may be written as 

e- . * 


£ dV l °ij d E ij + i j 7 \ dv = L dy / 

v» n 


e. . 
1 3 


cr • • d e . . 
13 


r — * t * t , u 

+ / 2 u i u i dV " /F - u. ) dV 

\T x 1 x 1 


(3.2.5) 


which finally reduces to 


a. . s 

/ dV / e ii da ii + / 7 \° u® dV >_ / T S u) dA 
¥ 0 if A x 1 


(3.2.6) 


Tlius for the material obeying the Drucker' s postulate of 
stability(3,4,5) , the principle may be written as 


f A 4 u i dA - f v J “i “i dy + / s ij( 0 i/)dV (3 - 2 ' 7) 


where !L = Any safe set of surface tractions 


/ 

¥ 


xj ij 


dV 


= Complementary energy associated with 

g 

stress field a . . , which does not 

(3 

violate yield condition. 
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a) SIMPLY SUPBDRTED BE M ON ELASTIC FOUNDATION : 

A unifom simply supported beam of elastic material 
resting on a Winkler model of foundation hawing a foundation 
constant K, is subjected to impact loading such that initial 
velocity distribution is known as shown in fig. (3.1b). The moment 
curvature relationship of the beam is as shown in fig. ( 3 . 1 c), 
treating the material elastic and perfectly plastic which also 
satisfies the stability postulate [ 3 , 4,5 ]. 

If the central transverse displacement of the dynamically 
loaded beam [figure (3 • 1 a,b)j at any time during the impulse is given 
by 5 , then to satisfy Eq,.(3.2.7) 

P.« £ K° + C S (5-2.8) 

neglecting all other terms including shear deformation. 

Where 

K° = Initial Energy terms (body forces not included) 

for initial velocity distribution of v(x) 

L /2 

= 2 j •§-. P A. v(x). v(x) dx (3.2.9) 

o 

which is analogus to term / -g-. p u.° . u.°. d¥ 

Y ix 

in Eq..( 3 «2 .7 )• 
s 

C = Complementary strain energy stored in the beam and 

foundation 

P represent the choice of T? as shown in Fig. (3.1a). 
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So C‘ 


_1_ , L , dV 2 . , i , 

2 El J ' 2^ + 2 ^ j Vv dx 

/-N ^ * 


dx 


(3.2.10) 


which is analogies to I \ .(o ± . S ~Jd¥ in E q. (3.2.7). 


ho W considering any specific value of P, the optimum 
hound con he obtained by maximizing <5 , such that 


35 

sF 


= 0 


(3.2.11) 


For this simply supported beam 0 < x <L 
admissible state, with a central load P, 


, in a statically 


Bending Moment = M 


x 


P_ 

4x 


X x 


(3-2.12) 


Deflection = W = ~ A 

2L Ax 


(3.2.13) 


Maximum Bending Moment 
at centre 

where, 


P 

4X * * 


(3.2.14) 


c. 


1 


Xx Coda X Lf. Cos XL 


A 


1 


Xx 


Cosh X LfCosX 1 


"f Cos X x gin X (Zr-x) - 
Sinh Xx Cosx(l-x) - 
Sin X x Co sbX (L-x) + 

Cos pc Sinh X (l—x) } (3.2.15) 

{ Cos he Sinipi (L-x)- 
CoshXx Sinx(L-x)+ 

Sin xx CosX(L-x)- 

SinhX x Cosx(L-x) > (3.2. 16 a) 
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rf> = Sinh A L + Sin A I, 
Cosh XL + Cos X L 

1/4 

in which X = [K/4EI ] 

So strain energy stored in the beam 
and foundation 


where 4 . = s 
1 X 


1 _ 

X 2 


° 128 EE ^ V * 1 ) 

l A* ^ d * 

t A,/ dx 


(5.2. 16b) 

(5.2.17) 

(3.2.18) 

(3.2.19) 

(3.2.20) ' 

energy can he more 


However, this complementary strain 


easily computed considering the desired solution Eq.(3.2.13) 
in the form of trigonometric series for the beam of finite 

length. Hateyai [9 ] has illustrated that a Sine series of the 
form 



2 Sin ftirx 
'h ° L 


(3.2.21) 


will satisfy completely the boundary condition at both ends. 


Strain energy of bending of -fee beam 


V = C s 
1 BEiM 


4 

El 

4L 3 


I 

n=1 


4 2 

n a 
n 


(3.2.22) 


l.l.T. f'/'IVFUR 
yEN T R A {_ LIBRARY 

^cc. No. /• • 
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and that of foundation 


V = C 

2 Foundation 


KL 


7 I a 2 
4 n.1 n 


For srna li v ax i at i on 


in 


a n ’ ^e re suiting change in strain 


will he 
9V 


energy 


av, 


. da = 

Hi n H 


3V 2 

+ r ,i. 


d2 » ’ ~ ,ia n 


n 
4 

_ TT El 4 FT 

- — n a dsu + £\r a ^ n 
2L 3 ^ °n 2 n H 


nu the sane time, because of the change da , the 

n 

loading P, a t a distance L/2, from the 


work done by 


support, is given by 


‘ v d a_ = P Sin d 


n 


nuc 

n a 
L n 


(5.2.23) 


Equating this change in potential energy to the change in the 
strain energy of the system 


4d- 


P sin n7rc - 1 EI 4 

v sm — - — = — — n 


0 , KL 

r •7 li S- + "T mmm 3. 

t' 9 t 3 n 2 


n 


So ’a^ 1 becomes, 


where c = L/ 2 . 


a = 
n 


2 PL^ 

tt 4 ei 


Sin 


n 7 r 


v, 4 KL'~ 

n + — 

tt nl 


Hence 


^ s 
G « 


128 El 




(3.2.24) 
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where 


“ 128 4 , Sm n ^/2 _ 2 

' 1 = l — T n L C -T ~ X 1 

n=l it 4 n 4 +4(AL/7r)' 


(3-2.25) 


“128 43 mr/2 2 

'1 = l ~J CAL) L [ 


n=l n 


7 

A ti 


n * + 4 ( XL/ tt) 


(3.2.26) 


and 


■?dt^ “ Sin m rr /2 Sin mrx/L 

W = ^ l [ — 

rEl n=l n 4 + 4(A L/rr) 4 


(3.2.27) 


Considering a state when "bending moment anywhere in the "beam 


is loss than the limiting moment M , 


P 

4 A 


< M. 


(3.2.28) 


as long as Ecj.. ( 3 • 2 .28 ) is satisfied, the entire "beam is in elastic 
range. How optimizing ' 6 ' with respect to P, in Eq,.(3.2.a), 

|r-(P. 6 ) = -L(K° + C S ) 


3 P + 6 - a p 


l£ s k°+ c s 

3P‘ P 


(3.2.29) 


s 

patting the value of C from Eq.(5.2.18) 


p = ( 128EI K^ 
<SM + M 


(3.2.30) 
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Ihus value of p, is given by Eq. £ 3 . 2 . 30 ) as long as H <K 

X“ 0 

4 X M 

_ o 

* ( 3 . 2 . 31 ) 


So P < 


i.e„ 


^128EI A 2 4XM 
l T~T , — ) < 


(3.2.32) 


So when tic response is elastic for the entire beam , the take of 

P given by Eq. ( 3 . 2 . 30 ) ciicmld be adopted and othenri.se when 

E<1-(5.2.31) is not satisfied, the valne of P given by the equality 

id Ed- (3.2.26) should be taken. Thus when the entire response 
I s e la stic 


5 < 


A C S 


« £ I 


On the otherhand if 


s°(VV* 


32 El 


■] for M > [ 
**' 


8 El E° <f > 2 i 


x (<A ^) 


3 (3.2.33) 


M < r 8Si k° 4* 2 ,s 

° ~ .2/. , a J 

A ( ^) 

Talcing the value of P from (3.2.28), 


4 X IJIq 


K° <P M 0 X (i>1+ ^ ) 


4 X M 


32 El .(j, 


Then 6 < 


+ 


(3.2.34) 
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Eq.(3.2o3) and Eq, ( 3 . 2 . 34 ) gives the general solution for a given 
value of initial kinetic energy input. As an illustration for 


comparison of results taking a half f 


sine curve for initial velocity 


dx o uci oution [Ref . 19 3 > given by 


v(x) = v Sin 4L? 

0 L 


(5.2.35) 


If _ / a 2 2 TT X 

o 2 T o 

?;here m = PA = mass of beam per unit length. 
_nt._cuu.cing a dimensionless impulse parameter 


mEI v 

S = ^ (b/AL) 

2M 


0 / w._ n 2 1} 




(3.2.36) 


and a dimensionless displacement parameter ' 


D = 


M q .6 


? 2 
mL ' v 
o 


(3.2.37) 


rni 


hen eq 


equation ( 3 . 2 . 33 ) reduces to 


M 6 
0 

T 2 2 

mL y 

0 


16 (s)* 3 ( ^ for s <_ 


and equation ( 3 * 2 . 34 ) reduces to 


(3.2.38) 


M .6 

0 - / 1 


mL v 


o < / 1 1 N v 

2 2 - ( 16 + 64S } Tl forS -4 


(3.2.39) 
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Hie values of S,D, for XL = 1.0 is given in Table ( 3 . 1 ). And S 
and D have been plotted on a semi-log scale for different values 
of XL in fig. ( 3 . 2 ). 

3.3 FEES -FREE BEAM ON ELASTIC FOUNDATION? 

The beam 0 <x< L with load P fig. ( 3 . I'd) in statically 
admissible state, has the following characteristics 


B.M. at any distance x = M = — C, 

x 4A 1 x 


v\ 

Deflection = W = A 

2*- AX 


(3.3.1) 

(3.3.2) 


Maximum B.M. at centre 

Maximum Deflection 
at the Centre 

H where 

c = _J, 

Ax SinhxLt- Sin XL 


= (P/4X)$ (3.3.3) 

= ^ C°sh>L + CosxL +2 _ n 

2K SinhxL + SinXL ' UO * 4; 

{Sinh A x [Sin Ax - Sin a(L-x)]- 
Cosh Ax[Gos Ax + Cos }(L-x)3+ 
Sin Ax [Sinh Ax - Sinh a(L-x) ] 
Cos Ax [Cosh Ax+ Cosh a(L-x)]} 

(3.3.5) 


A Xx Sinh XL + Sin XL { Cosh Xx CosX ( L- x) + 

Cos Ax Cosh X (L~x) - 
Sinh Ax Sin a(L-x) + 

Sin Ax Sinh a(L-x) + 

2 Cosh Ax Cos Ax } (3.3*6) 
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Cosh XL - Gos-AL 
Sinh lL + Sin \L 


( 5 . 5 . 7 ) 


So with free-free boundary conditions, the only change from simply 

s 

supported Beam Analysis, will he "that the C will he given hy 



128EI 




(3-5.8) 


whore 



/ (0 x / dx 


( 3 . 3 . 9 ) 





dx 


(3.3.10) 


Aid C and are given by Eq.( 3 .3.5 ) and (3.3 .6). 


Following the analysis presented in the previous section 
(Sec. 3.2), the final equations governing the deflection can he 
obtained as 


M .5 1 

— 2“ 2 < (1/S) 2 (<^>l) for S <i (3.3.11) 

ml v 

0 


M . <$ 
0 

_2 2 
ml v 
o 




JL 
X 1 


for S > i 


(3.3.12) 


She value of S and D for Xl =1.0 is given in Table (3. 2). 


S and 1 have been plotted on a semi -log scale for different value 

of i* 1 I'ig. (3.3). The values of 1 for free-free case are 

much less than for simply supported case. 
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54 ILHJSTEMED EXAMPI E % 

(a) As an example of the application of this principle 
of displacements computation, a simply supported Beam on Elastic 
foundation subjected to impulse loading is considered and the 
computed bounds are compared with the exact solution obtained 
earlier (Sec. 2.l) using Modal Analysis. The accuracy of computed 
bounds is reflected in this comparision. But for a very simple 
case with -well defined impulse, exact solution are available and 
so this comparision does not necessarily reflect the accuracy of 
these bounds for different boundary and initial conditions. 

Considering a uniform simply supported beam on elastic 
foundation -with the. initial condition, 


W( x, 0) = 0 


0¥ 
3 t 


(x, 0) 


Sin tt x 

Y 

o L 


(344) 


Hie governing differential equation is 

EI lihid + a jkad + ra(x>t ) , 0 

3 x 3 t 

has the solution, 

(A Cos w t + B Sin w t) 
v n n n n ' 


v7(x,t) = l 
n»1 


Sin 


n tt x 


(34.2) 
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Using the initial conditions given by Eqns.( 5*4 # 1 ) 


v 

¥(y.,t) = ~ Si n If Sin 


(54 .5) 


Thus showing that the bean will vibrate only in first mode 
where the first natural frequency is given by 


4 ^ 4 "^ 

L El TT 


m 


y 1/4 

X ~ t AEI^ 


Eq. (54.3) becomes 


¥(x,t) = 


v L' 
0 


7T 


o • tt^- 

.2 ° in ~L 


Sin w^t 


/1+ 4l^L/vy 


[— ] 

L EI J 


1/2 


For maxLnum central deflection 


Sin 


•n-x 


= 1.0 


= 1/2 


Allowing tills to happen at time t* , such that 


Sin w^t' = 1 


W 


v o L ' 


max 


_ — i — • 

TT . i. A { 1 T L 4 -1/-2 


1/2 


[ 1+ 4 ( 1 V 1T 3 3 " 


For XI » ir . 



n 1 /' 

if . \ = 0.(4550 v l . r =-] 

(.maximum) o L El J 


(34.4) 


Using Eq,.(3.1.33) to evaluate the hound for central deflection. 


0 = 


t K (<h,+ V, 

^ 32EI 3 


(34 .5) 


where 


KT 


mlvc 

4 


d>^ and are obtained from Eqs. (3.1.25) and ( 3 . 1.26) for XL air, 
tricing first three terms only, as the series is rapidly 

converging and are obtained as 


^ + *, - 0.2572 If 


(3*4 . 6 ) 


Hence from Eqns. (34.5) and ( 3.4 ; 6 ) 


5 = v 0 .L 2 . O.O 515 [^} lU (34.5) 

Comparing this with exact solution given by Eon. ( 3 . 4 . 4 ) the computed 
bound is only 12 $ higher than the exact value. 

(b) Considering simply supported beam, resting on an elastic 
foundation, which is subjected to a central impact with a mans 
and velocity v . Assuming no loss of energy, the kinetic energy 


incut is taken as 
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Using Sqn „(3.2.33), upper "bound on displacement 


x s 


max 


K ( vm+ % ) i 

[ 1 J_i 

32EI J 


Hence 




m 2 <‘ 1 + ^ 2 


max. 


32 El 


] • v 
J n 


1 

8 V o 


V^V t) 

[— 1-1 

*■ El J 


Tricing w ^ = 0.2572 L 5 


5 = — t [ 

max. 8 o l 


1 m x 0.2572 ^ 3/2 

r _ Z v 


El 


] *L 


F °r v o = 100 inch/sec., m 2 - 3.78 , El = 6 x 10 5 lb-inch 2 , L = 40 


o — 0.22 6 = 0.101 inch 

max ‘ 1 xnL,n * 


Thd maximum displacement by modal analysis 


Vx “ °- 077 i”<* 


Thus the computed bound is 1.30 times more than tho 
maximum value given by Modal Analysis. 
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TAB IE 3.1 


FOR SIMPLY SDPIOECED BEIM 
ON ELiSTIC 
FOUNDATION 


XL 

4 

4>/AL 

0.5 

0.60 

1.20 

1.0 

1.00 

1.00 

2.0 

1.35 

0.675 

IT 

1 .09 

0.348 

4.0 

O.99 

0.248 

5.0 

0.982 

O.196 

2tt 

1.001 

0.158 

8.0 

1.01 

0.125 

For 

= 1.0, Value of S and D 

s 

I) 

: REMARKS 

0.10 

0.187 

ELASTIC 



HINGE 

0.20 

0.140 


0.25 

0.125 



— - - - ~ 


0.50 

0.0938 


1.00 

0.0/82 


2.00 

O.O7C4 


5.00 

0.0657 


10.00 

0.0641 


100.00 

• 

0.06 27 


• 

• 

CO 

0.0626 




TABLE 3.2 
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FOR FREE FREE BE AM OR ELASTIC FOUHDATIOF 


XL 

* 

$/XL 

0.50 

0.252 

O.5C4 

1.00 

O.5O5 

0.^05 

2.00 

O.924 

O.464 

ir 

I.O9 

0.547 

2rr 

O.995 

0.160 

8.0 

0.98 

0.124 


For X L = 1.0 

, Value of S and D 

s ! 

D 

REMARKS 

0.10 

O.O99 

ELASTIC 

0.20 

O.O7O 

RANGE 

0.25 

_ _0. C62 


0.50 

O.O47O 


1 .00 

O.O39O 


2.00 

O.O35O 


5.00 

O.O325 


10.00 

O.O32O 


100.00 

O.O312 


• 

* 

CO 

0.0310 
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3.5. BOUNDS ON DISPLACEMENT AND ENERGY USING VARIATIONAL PRINCIPLE. 

An attempt is made to evaluate the ipper bound on energy and 
displacement for beams on elastic foundations subjected to a 
general class of dynamic loading q(x,t). The governing differential 
equation of motion is (neglecting the damping) 

2 2 2 

— 2 [ E I Cx)-~|] + m(x) —•+ K(x) W = q(x,t) (3.5.1) 

Defining EI(x) = S(x) = Bending Stiffness. 

With initial conditions as 

3W 

W(0,t) =0 and — (x,0) = 0 (3.5.2) 

The following end conditions of the beam are 

(a) Pinned - Pinned 

W(0,t) = 0 , S(0) — = 0 

3x 2 

(3.5.3) 

W(L,t) = 0, S(L) — , o 

3x 2 

(b) Clamped - Clamped 


w (0,t) = 0 , || (0,t) = 0 

W(L,t) = 0 * H CL,t) = o 


(3.5.4) 



65 


4 

(c) Clamped - Free 


W(0,t)= 0 , 


3W(0,t) 

3X 


0 


C3.5.5) 


S(L) „ 0 
3x 2 


ax 


CSCL) = 

Im 

ax 


(d.) Free- Free 


S(0) 


3 W(0,t) _ Q 


3x z 


S(L) 


3 2 W(L,t) 

ax 2 


a_ 

ax 


1_ 

ax 


[SCO) 


[S (L) 


a 2 w(o,t) _ 

ax 2 


a 2 w(L,t) = Q 
ax 2 


(3.5.6) 


Considering an elastic beam resting on elastic foundation (Winkler 
Model) with one of the above mentioned boundary conditions, * x' 
is the axial distance along the beam, x = 0 , being the one end and 
x = L the other end, * t* is the time parameter and VJ(x, t) is the 
lateral displacement. The bending stifftiess is denoted by EI(x) = S(x) 
and the mass per unit length by m(x) . S(x), n(x) assumed to be 
symmetric about x = L/2. 

If the dynamic load q(x,t) is applied to the beam at t = 0, 
the equation of motion is given by 

iL [EICx) + mC x) ilS&lL + KW(x,t) = q(x,t) (3.5.7) 

3x Z 3x Z 3x 
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Neglecting shear deformation and rotatory inertia. defining the 
function 


V* = 2 [Strain Energy + Kinetic Energy] 


(3.5.8) 


r* = £ 

0 


[S(x) 


f 3 2 W 2 
3x z 


+ KW 2 + m(x) 


(— ) 2 ] 
3t' J 


dx 


(3.5.9) 


The time rate of change of V*during motion of the beam is given by 


dV* 0 f L ro , , 3 2 W 3 3 W 

j— = 2 / [S(x) — ~ T ~ 
at 0 3x 3x Z 3t 


. m 3W ^ . . 3W 3 2 W. , 

+ KH + ^ dx 


Integrating by parts, together with boundary condition, the above 
equation can be expressed as 


dV* 

dt 


L 2 2 \ 7 

= 2 J [£-. (S (x) ■“■)] + KW * m(x) 
0 3x Z 3x 


3 2 W. 3W 
3t 2 3t 



With the use of equation of motion given by Eq. (3.5.7) 


dV* 

dt 


= 2 / q(x,t) ~ dx 
0 3t 


(3.5.10) 


It follows from Schwarz inequality [11] and definition of V* (Eq. 
3.5.9) that 


dV* 

dt 


< 2 a(t) 


L 

{/ m(x) 
0 


A 2 


dx ) 


1/2 


<2 a(t) V* 2 



Defining V* = where Vj >_ 0 


dV* r dV l 
dt 1 dt 


/ <2 aft) . V_ 

dt ” 1 


< a(t) 


(3.5.10a) 


Integrating Eq. (3.5.10a) and using the initial conditions Eq. (3.5.2) 


Vj / a(x) d t 
0 


V* < / a(x) dx 


(3.5.11) 




{ f ( S ~~p 2 + ra Cx) c ~) 2 + KW 2 ] dxF < / a(x) dx(3. 5. 12) 

C 3x - B 0 


Tfiis inequality furnishes an upper bound on energy. 

For obtaining a bound on displacement, it is sufficient to derive an 
inequality of the fom [26,27] 


/ L r B \2 . 2 .,2 

J (—j) dx > E> ft 
0 d x l — m 


(3.5.13) 


Inhere W = Maximum Iwl for 0 < x < L 


and B = constant 
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or 


/ (W xx ) 2 dx > B 2 W 
© 


Equation (3.5.13) is to be satisfied at time *t’ by all possible, 
displacements W. In addition the largest value of constant B in 
Eq. (3.5,13) is to be obtained. 


The value of f B* is derived by the procedure suggested by 
Plant [26], As stated above, the above inequality (Eq. 3.5.13) is to 
hold for any time t, the problem of determining the extremum value 
of constant B; can be reduced to a time-independant variational 
problem [26,27], 


To solve the variational problem, an auxiliary function u(x) 
is introduced. The function u(x) will be called admissible if it 
satisfies the appropriate boundary conditions for the beam, resting 
on elastic foundation and if u(x) is continuous and if 
L 2 

/ (u ) dx < 00 (3.5.14) 

0 

Considering all admissible deflections for which 
L , 7 

/ (u vv r dx = r (3.5.15) 

0 x 

relation given by Eq. (3.5.15) may be viewed as a normalisation 
condition on admissible deflections. The largest value of ’B' is 
then determined in two steps. First, the particular admissible 
function which satisfies Eq. (3.5.15) and which yields the maximum 
deflection at some arbitrary point along a beam on elastic foundation. 
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say x - a is found. Then, the particular value of 'a* say a = a* 
at which the foregoing set of maxima has an absolute maximum, is 
determined. Once the location a* (therefore the maximum deflection), 
is known, the extremum value of B, may be calculated. 

In mathematical terms, the first step is to maximize the 
function u(a) for any point 0 a _< 1 , subject to subsidiary 
condition given by Eq, (3.5,15). 

The problem is clearly suited for application of Lagrange’s 
method of undetermined multipliers. The Lagrange’s function can be 
written as 

L 

F(u) - u(a) + X* {/ (u xx ) 2 dx - C 2 } (3.5.16) 

In Eq. (3.5.16), u has to be determined such that F(u)is stationary, 
where X* is a undetermined multiplier. 

Using these, the final expression for u(x) can be written 
[26,27] as 


"W = ~p~z g(x[a) 
^(aja) 


(3.5.17) 


Where g(x| a) , is the influence function for the beam, resting on 
elastic foundation and X* has been determined in terras of C from 
Eq. (3.5.15). 
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influence for the bean, resting on the elastic 

foundation is defined in the custonary fashion, as the displacenent 
at *x' due to a unit lead at 'a'. From Eq. (3.S.17) it can be 
concluded that for a given strain energy, the maximua deflection at 

* = a “ ^ved by applying a concentrated force of appropriate 
magnitude at that point* 

The second step is to determine the value of a* at which the 
magnitude of u(a) and hence g(a|a) is maximum. 


The point a* is either an end point 
which satisfies 


or an interior point 


dg(a|a) _ 3g(a|a) a £ ( x |a) 

= [ + ] _ Q 

da 3x 3 a x=a 


(3.5.18) 


Using Maxwell's Reciprocal theorem g(x|a) . gta|x), Eq. ( 3 . 5 . 
also be written as 


18) can 


3g(x|a) 

"^x=a = 0 


(3.5.19) 


Which implies that u x is zero at the point of application of the 

After solving for a*, the maximum deflection at any point 
can be written as 


u(x) 


•^(a^a*) 


g(xj a*) 


(3.5.20) 


and the absolute maximum itself 


as 



71 


u m = c ^gC a *| a *) (3.5.21) 

The extremun value o£ B is then given by 

/ ^xx ^ 2 dx 

2 C xx 1 

B z = = (3.5.22) 

u g(a*|a*) 

The quantity g(a*|a*) is the maximum possible displacement of the 
beam, resting on elastic foundation under a unit concentrated load. 
For 


(a) Kinged-Hinged ends 

a* = L/2 

Sinh XL - Sin XL 

g(a*|a*) = ~r (3.5.23) 

Cosh XL + Cos XL 

" heie - x= [ 4ir ll/4 - 

(b) Fixed- Fixed ends 

a* = L/2 

. Cosh XL + Cos XL-2 

g(a*| a*) = w (3.5.24) 

' v Sinh XL + Sin XL 

(c) Fixed- Free ends 

a* = L 

. Sinh2XL - Sin 2XL 

g(a*|a*) = f — 

Cosh z XL + Cos^ XL 


(3.5.25) 
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(d) Free-Free ends 


a* = L/2 


g Ca* | a*) = 


^ Cosh XL + Cos XL +2 
2K Sinh XL + Sin XL 


or g(a* | a*) 


2x Cosh XL/2 Cos XL/2 
K Sinh XL + Sin XL 


, depending upon the value of XL 


Putting g(a*|a*) = u g , bound on displacement can be obtained as 


/ S(x) C^|) 2 dx > W 2 

o ax z u s 


(3.5.26) 


Combining Ea. (3.5*12) and (3*5*27), the bound on the displacement 
can be written as 

1 

7 t 

|w(x,t) I £ Ug / a(T) dx for 0 ■< x <_ L (3.5.27) 

0 

2 A 

L q (x,t) 7 

where a(t) = [/ * — dx]“ 

0 n(x) 

If loading can be separated in the form 


q(x,t) = <{>(t) f(x) 

Eq. (3.5. 27) can be written as 

1 

y t 

iw(x,t)|_<g s / U(T)j dx 0 < x < L 
0 “ 


(3.5.28) 


(3.5.29) 
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where 


8 „ 


= u s / 


0 


^(x) 

ra(x) 


dx 


Thus for general load q(x,t) Eq. (3.5.27) furnishes the 
x^sper bornd on displacement and for the load q(x,t) of the form 
given by Eq. (3.5.28) the inequality given by Eq. (3.5.29) provides 
an upper bound on the response of the beam resting on elastic 
foundation. 

3.6 IMPACT LOAD. 

If the q(x,t) is given by Hertz Law of impact, then 


q(x,t) = m 2 v 0 (l+e) —■ Sin ~ for t < ^ 


= 0 


for t > T t 


( 3 . 6 . 1 ) 


For free-free beam on elastic foundation. 


— , Cosh XL + Cos XL + 2 . 

ul = [± .] 1/2 

S toy J 

Sinh XL + Sin XL 


Then bound on displacement is given by 


, Cosh XL + Cos XL + 2 , n t 

|W(x,t)| = ” U 


Sinh XL + Sin XL 


■ir' / a(T) dr 
0 


Putting the value of q(x,t) from Eq. (3.5.1) 


a(r) = m 2 v Q (1+e) Sin ^ 
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t . m 2 v Cl+e)ir T t 

SO / o(T) dr - JhhK / L 

✓V O r P » 


— / L Sin J£ dt 

2 Tt 0 T l 


a Cosh XL + Cos XL+ 2 i /? mo V-fl+e) 

Il»tx,t)| = fM -] 1/2 — 


Sinh XL + Sin XL 


/L/pA (3.6.2) 


For a beam with the following parameters 


X = ri.080, *L = 3.20, L = 40 inch , m 2 = 3.75 

v q = 9,97 inch/sec., e = 0,220, p = 0.29 lb/inch^ A = 1.0 inch^ 

So |VJ(x, t) | = 29. 8 * 10 J inches-. 


This is about two times as predicted by Modal analysis and 
three times as predicated by Energy Method. 
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3.7 RESULTS AND DISCUSSION 

The displacement bound within the elastic range is given by 
Eq. (3.2.33) and the permanent deformation is given by Eq. (3.2.34). 

For a half sine impulse, the displacement parameter, D and impulse 
parameter S, have been shown in Figs. (3.2) and (3.3) for simply 
supported and free-free beam on elastic foundation for different 
values of XL. These values are also given in Table (3.1) and (3.2). 

The Fig. (3.4) shews the impulse and displacement parameter for a 
Simply Supported beam, given by Martin [19], and Seiler, Cotter and Symonds 
[30], It is observed that for beams on elastic foundation, comparatively 
for a smaller value of impulse parameter S, the response is in elastic 
range as compared to beams without foundation. For a half sine 
inpulse, v(x) = v q Sin mx/L , the upper bound on displacement (within 
elastic range) given by Eq. (3.4.5) is 12% higher than the exact value, 
evaluated by Modal Analysis, as in Eq. (3.4.6). For X° = m 9 v^ , i.e. 
a central impact by a mass m 2 with a velocity v c the computed bound 
is 1.30 times the maximum displacement given by Modal Analysis for 
Simply Supported beam. 

The displacement bounds for a dynamic load q(x,t) using 
variational principles has been given in Eq. (3.5.28) and Eq. (3.5.30). 
These boundshave been computed for Hertz impact force given by 
Eq. (3.5.2). For a free-free beam, this is about two times more 
than that predicted by Modal Analysis and three times as predicted 
by Energy Method, 
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CHAPTER 4 


PROPAGATION OF WASTES IN BEAMS ON ELASTIC FOUNDATION AFTER 

TRANSVERSE IMPACT 

4.1 INTRODUCTION 

The study of wave motion of a beam on elastic foundation 

0 

undergoing flexural vibration after impact is presented here to 
visualise how the flexural pulse and energy propagate in the beam. 

The medium can be classified as dispersive or non-dispersive depending 
on the propagation of waves. 

4.2 BASIC EQUATIONS 

If the shear deformation and rotatory inertia effects are 
neglected, the differential equation for the free transverse 
vibration of the uniform Euler-Bemoulli beam on Winkler Foundation 
is given by 


3x 4 


PA 

El 


3 2 W(x, t) 

at 2 


EI 


w(x,t3 = o 


( 4 . 2 . 1 ) 


Assuming that the solution of Eg.. (4»2,l) is a simple 

» 

harmonic wave travelling with velocity v in the positive in- 
direction, so that its form is 


W(x,t) = B' Cos -§~ (x-vt) 
0 


( 4 . 2 . 2 ) 


Where 8 is the wave length and B' is a constant. Substituting 

Eq. (4.2.2) in (4.2.1 ), the velocity of propagation of Sinusoidal wave 
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v 


/ 51 / ! + _ol 

e pa ei • 16 „4 


(4.2.3) 


Miere v = phase velocity. 

For K = 0, i.e. for beam without foundation 


v 


2tt / EI 
0 PA 


(4.2.3a) 


So the waves move faster in beams on elastic foundation. Therefore 
the velocity of propagation of Sinusoidal flexural wave is not- 
constant but varies with wavelength as given by Eq.. (4.2.3)* Usually 
the velocity of propagation v of a Simple Harmonic wave is Called 
wave velocity or phase velocity. 


Anon-harmonic flexural, pulse £23, 24 3 maybe regarded as 
consisting of a super position of harmonic wares of different 
wavelengths. Each of these waves have a different phase velocity, so 
it follows that a fluexnral wave of arbitrary shape cannot propagate 
alcng the beam without dispersion, which resilts in a change in the 
shape of the pulse. A medium exhibiting a wave velocity depending 
upon the wavelength is called a dispersive medium. Hence, this 
uniform beam on elastic foundation is a dispersive medium for 
propagation of flexural pulse. The only flexural, wave profile 
propagating in a uniform beam wilhout alteringits shape is the 
simple harmonic wave. A pulse consisting of a group of harmonic uasre s 
is called a wave packet [23] . ind the velocity with which such a 
group of waves is propagated is called Group Velocity, denoted by U. 
The group velocity is the velocity with which the energy is propagated. 
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Considering a we packet consisting of two simple harmonic washes of 
equal amplitudes, of wavelengths 6 and e + A 6 and of wave velocities 
v and v + Av, respectively, the time interval necessary for crest «e‘ 
to catch up with crest E, as shovm in Figure (4. l) is denoted by t. 

The new common positions are denoted by e f and E 1 , respectively, lave I 
travels with velocity v, wave II with velocity v + dv, and the packet 
with group velocity U [34] . From the Figure (4.l), 


A6 + (v + Av) T = VT 

0 + VT = Ur 


(4 «2«4 ) 


Combining Eq. (4.2.4) and letting the wavelength difference become 
very small, the group velocity expression can be obt ain ed as 


U = v 



(4.2.5) 


The Simple harmonic wave described by Cos — (x-vt) is called 

0 

the carrier wave and the envelope shown in Fig. (4.1) represented 
by 2B* Co s (~) (x-Ct), is regarded as moving with group velocity IJ. 


Again considering the propagation of flexural waves and 
substituting Eq. (4.2.2) in Eq. (4.2.5) 


U 


v - 0 


dv 

d0 


= v {1 + 


A 2 El 
S } PA 


JL (V 

PA 12 it 


f 2 tr 2 El K . 6 .2 
^8 ' PA pA ^2 tt 


For K » 0, U = v (1 + l) = 2v. 


( 4 . 2 . 6 ) 
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But for beams on elastic foundation obviously 


U < 2v (4.2.7) 

As the wavelength 9 approaches zero, both the wave and group velocities 
approach infinity. The group velocity is the velocity with which 
energy propagates. So it can be concluded that for a pulse consisting 
primarily of infinitely short waves, the energy is propagated with 
infinite velocity, which does not seem physically possible [J4] * 

The reason for this puzzling result can possibly be explained 
by the fact that for very short waves one cannot neglect rotatory 
inertia. When wavelength is of the same order of magnitude as the 
depth of the beam, the rotatory inertia effect is of the same order 
of magnitude as the translational inertia effects [34] . 

So neglecting shear deformation but considering the rotatory 
inertia effects in the formation, the differential equation of motion 
for a uniform beam on elastic foundation becomes 


3 4 W pA S 2 W 
3X 4 E * 3t 2 


PA 


a 2 w 
2 2 

3x 3t z 



0 


Ihere K 
S 


radius of gyration. 


( 4 . 2 . 8 ) 


Assuming a solution in the form given by Eq.. (4.2.2) the wave velocity 
can be obtained as 


v = ^7 p“ {1 + 


4tt 2 K 


1 

2 


n K i e 1 
+ EI i6* 4 


(4.2.9) 


111 ere rw is the wave propagational velocity for the longitudinal 
impacts 



84 


For K = 0, i.e. for beams without foundation 


v= /e/p {1 + — | 1 -1 ^ 2 

Similarly the group velocity becomes 


(4.2.10) 


u = Se/p~ {1 + 


e 2 -1/2 k 6 4 

} X {1 + — — } X 

4ir 2 K ' EI 16r 4 


Q 4 - 3/2 

1 - 4(1 + — : V' 

417 K g 2 

[ + 6 ( 1 


, , EI 16 ir 4 "s 

(1 + — — — / 

K 0 4 


6 2 + 4tt 2 K 2 

o 


-)] 


For K = 0, 


(4.2.11) 


u = sm {i + — r— : r 1/2 a + — - 

c g 


4ir 2 K tj.,^ 2 v 2 fa 2 


l+4ir* K“/0 

g 


( 4 . 2 . 12 ) 


Hence in Eq.. (4.2.1 1 ) and ( 4 . 2. 12), it can be seen that when 6+0, 
group and wave velocity tend to the constant value / E/p . 


4.3 EESULTS JTD DISCUSSIOH 

The v/ave and group velocity for beams on elastic foundation 
and for beam without foundation has been plotted in Fig. (4.2). For 
beams on elastic foundation the flexural wave velocity is also less 
than the group velocity. The group velocity with which eneigy propagates 
for beams without foundation is twice the flexural wave velocity, 

•whereas for beams on elastic foundation, it is less thai twice the 
flexural wave velocity. However these differences are only 
pronounced for higher wave lengths. Theflexural wave moves faster 
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in a beam on elastic foundation than in abeam without elastic 
foundation, whereas energy moves at a slower rate along a beam on 
elastic foundation than along a beam without foundation, which is 
shown in Fig. (4.2). The beam on elastic foundation is a dispersive 
medium for wave propagation, as wave velocity is dependent on wave 
length, as in Eq. (4.2.3). 
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TABLE 4.1 



& 

El ~ 

0.10 , El/ A = 

1,0 


WAVE 

LENGTH 

WAVE 
VELOCITY 
FOR BEJM 

WAVE VELOCITY 
FOR BEAM ON 
ELASTIC FOUN- 
DATION 

GROUP 
VELOCITY 
FOR BEAM 

GROUP 
VELOCITY 
FOR BEJM 

ON ELASTIC 
FOUNDATION 

1.0 

6.280 

6.281 

12.56 

12,5000 

2.0 

5.140 

3.142 

6.28 

6.2500 

4.0 

1.570 

1.573 

3.14 

3.1000 

5.0 

1.255 

1.256 

2.510 

2.5000 

10.0 

0.628 

0.803 

1.256 

0.7720 

20.0 

0.314 

O.768 

0.628 

O.O568 
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CHAPTER 5 

EXPEEDEUTA! SET-UP M D ME ASU BEMEU TS 

5.1 IUTBODUCTION 

To study the response of "beam on elastic foundation, 
experimentally, a suitable set-up is needed to do so. The set-up 
in the Dynamics and vibration laboratory of Mechanical Engin ee ring 
Department, I.I.T. Kanpur is used for the experimental studies 
reported here. Hie details of the set up are shown in fig. ( 5 . 1,5,2) 
and in plate ( 2 , 3 , 4 , 5 * 9 ). 

5 .2 BRIEF DESCBIPTICN OP SETUJP 3 

in experiment has been designed to verify the theoretical 
predictions! 12, 13 , 32 , 35 ] The experimental arrangement is given in 
figure (5*1 * 5-2). Details of the set-up are shown in plates ( 2 , 3 , 4 , 
5,9 ). Ike hammer is supported by means of low stiffness 

stainless steel tube to minimize the energy absorption and also to 
ensure the movement of the striking mass in a vertical plane. At 
the top, the suspension block (fig. 5.2) has been precisely made 
to avoid frictional farces. Sliding arrangement at the top makes 
possible the adjustment of the suspension point of tine hammer 
whenever desired. Scale attached to the top supporting block 
measures (of the order O. 5 0 ) the angle of swing for the hammer. 
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Neglecting the friction losses and assuming pendulum to he a simple 
one, expression for striking velocity V Q is obtained by equating 
the maximum potential and kinetic energies' of the- hammer; and can 
he written as 

l. 

V q = 2gR (l- Cos 6) (5.2.1) 

where R = Length of Pendulum. 

5 . 3 EXPERIME mU, MOLEL: 

Seam of mild steel and of uniform cross-section was put 
on ten almost identical linear springs at a spacing of four inches. 
The details of the characteristics of the beams and springs used 
for the experimental study are given in Table 5.1* With a help of 
washer, nuts and bolts beam and springs were connected such that 
different teams can be attached to the springs. These springs were 
welded to a fairly heavy channel making the base rigid. Die ends 
of the beam were kept free-free. Hence a uniform free-free Euler- 
Bernoulli beam resting on Winkler model of foundation was simulated 
in the experimental set-up as shorn in Plate 1 . 

5.4 IN STHJMENTIiTION s 

In studying experimentally the response of the beam 
subjected to transverse impact, various parameters such as bending 
strain, deflection, acceleration, velocity of the beam have been 
measured. The very short duration of the impact (of ihe order 100 
micro seconds) needs effective design of the electronic circuitry 
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as well as the proper selection of the transducers for me asurin g 
various parameters involved. The details are discussed in the 
following sections. 

5.5 BERING CON STMT: 

All the ten springs which were used in the model were 
tested on a uni -axial compression testing machine. It was found 
that the load -Kef lection relationship was fairly linear during 
loading, unloading and reloading. After fabricating the model, 
it was pit on platform of a tri -axial testing machine (aIMIL). 

A .proving ring was used to measure the load at the centre of the 
beam and deflection, at centre (i.e. below the load pt), was 
measured with a dial gage (1.0 - O.Oi ■) . Hie deflection at 
centre under the load P, for a free-free beam on elastic founda- 
tion is given by, 

(5*5*1) 


(5.5*2) 


__ £ ^ Gosh + Cos 4- 2 

2K Sirh XL + Sin XL 


where 


or 


X* . 


K_ 

v^j+EI 

PL? ^ 
8EI(XL)’ 


Coda XL + Cos XL + 2 
Sinh XL + Sin XL 


The following observations were made 


P in lbs - W in mm 

50 - O.56 


100 


1.1 



92 


From Eq..(5 .5 »2), by trial and error method 

XL = 3.2 

So X = C.Cg inch 

K = 99 .8 Ib/in 2 

And width of "beam = 1.98 inches. 

Thus this model represented a foundation having a modulus - 5 0.0 Ib/in 
5 .6 3TBAI1T MEASTJEEMEHTS t 

Strain gages of hounded wire type have been employed for 
measurement of strain. At five places strain gages have be^n 
installed. At point of impact, two strain gages on the back side 
of the beam were installed and two dummy gages were installed, on 
similar material as that of beam . At other places one strain gage 
on each side of the beam were cemented. 

The leads from the strain gages were connected to Type 3C66 
carrier Amplifier (Bridge Balance) plug-in units. The plug-in unit is 
designed to be used with 564 -oscilloscopes. This unit is self 
contained and requires no external equipment other than the transducer 

For strain gage at the middle, four External Am Connection ?ras made 
as in fig. (5. 3). For strain gages at the other places, two External 
Arm Connection was made as shown in fig. ( 5 . 5 )* 

The type 3 C 66 plug-in unit was inserted into the proper 
plug-in compartment of the oscilloscope. Strain gages were connected 
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to the Type 3 C 66 . The Unit Input Bridge Circuit was Balanced, and 
with the strain gages in an unstressed condition, the bridge can 
be checked for proper balance by placing the FUNCTION switch in 
the BALANCE position. If the bridge is correctly balanced, no 
signal appears on the Cathode ray tube (CRT). 

5.7 STRAIN -GAGE CALIBRATION s 

The Type JC66 plug-in unit can be calibrated with the 

standard 15 OK calibration resistor supplied with the unit. This 

resistor is stable and accurate with good electrioal properties 

at 25 KC . Operating the Calibrate Switch on the type 3^66 unit 

from front panel which connects the calibration resistor in parallel 

with the bridge arms and EXT. ARMS Connector. Ibis reduces the 

resistance of the bridge arm, and is equivalent to subjecting the 

strain gage to a negative (compressive)strain And the equivalent 
strain is given by 

R x 10 6 

Equivalent strain = ( 5 . 7 * 1 ) 

< gj -)(h sg + W 


where 

G.E. 

- Gage Factor of strain gage 


^SG 

= Strain Gage Re si stance 


^CAL 

=t Calibration Resistor Value 

For 

G.F. 

= 2.05, R gG . 113.15 ohms, R Cil « 15CK 


Equivalent Strain = -37 0 . 
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Hie GillT Control is adjusted to make the S TEA TTJ /nTV 

% 

switch direct reading for strain measurements. 

In equation (5*7*l)» the strain simulated by the 15 Q£ 
calibration resistance for specific non-standard condition is -37O 
strain. Using this value, the type 3C 66 plug in unit was 
calibrated at 100 strain/div sensitivity. 

Division of Trace Deflection = "^05 = -7 divisions. 

The actual calibration was performed by 

1 . Setting th e STMIN/DIV switch to 1 00. 

2. Trace was positioned to top of CRT. 

3. Operating the 'Calibrate' switch, the trace 
is deflected downwards. 

4. GiilH Control was adjusted for a deflection of 
3.7O major gratitude divisions. 

Thus the STJLJN/DIV was calibrated for direct reading at all 
settings . 

5.8 DEFLECTION MEiSUEEMENTS : 

The deflection at various points were measured with the 
help of Vibration Pick up, (of piezo -electric ceramic accelerometer- 
type), Vibration meter and storage oscilloscope as shown in plate 


. The pick-up placement was accomplished with 156O-P35 
Permanent -Magnet clamp. This is a disk-shaped magnet with a 
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captive screw which was threaded into the tapped holes in the 
pick -Tip. With the FUNCTION switch at DISP, the vibration meter 
performs a double integration on the acceleration response of the 
pick-up. Further the deflection- time trace (plate 6 ) was 
obtained on the storage oscilloscope. 

Phe velocity and acceleration response of the beam was 
obtained with the FUNCTION switch of oscilloscope at VEL and 
ACCEL respectively. 

Hie calibration of the vertical scale of the oscillosoope 
was done with the help of an Electro -dynamic shaker and a Power 
Amplifier as shown in plate 

5 .9 SPECIFICATIONS OF INSTRUMENTS USED : 

1 * Oscilloscope 

a. 564 Storage Oscilloscope 

This is a special purpose oscilloscope designed to store 
Cathode-ray tube, and displays for viewing or photographing up to 
an hour after application of the input signal. This can he operated 
as a conventional oscilloscope also. There are two separate 
storage screens of the Cathode-ray tube. Either the upper or lower 
storage screen in storage mode can he used while the other can he 
operated in non -storage screen. 

Hie Plug-in units used are 

1* Type 3A3 Amplifier Unit. 

2. Type 3B4 Time Base Plug in Unit. 
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3. Type 3C66 Bridge Balance. 

TEKTRONIX INC . BEAVERTON, OREGON (U.S.A.) 

2 . Vibration Hck-DT 

v.' This is an acceleration type piezo -electric pick-up 
used to measure amplitudes with, the help of vibration meter. This 
can he held in position on any vibrating body to obtain the 
amplitudes, or it can be fixed to the vibrating body with magnetic 
attachment. 

GR Type - 15 1560 - P 52 

SERI IB NO 3397 

Pick-up sensitivity - 67*5 mV/g. 

Manufacturers: General Radio Company, Concord, Massachusetts, U.S.A. 

By Shore Brothers, Inc. Evanston, 111., U.S.A. 

3 . V ibration Meter 

The type 1553-AK Vibration Meter is a portable, gensre„l 
purpose instrument for measuring mechanical vibration when used 
with -Hie Type 1560-P52 Pick-up it is direct-reading in 

Acceleration ( inch/second ) 

Velocity (inch/second) 

Displacement ( milli inches) 

The vibration Meter consists of a piezo electric ceramic accelerometer 
-type pick-up and a sensitive vacuum-tube voltmeter. It has 
integrating circuits to convert acceleration to velocity or displacement 
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signals, a differentiating circuit to convert acceleration 
signals to jerk signals. 

Vibration Meter Serial No. 
1557-0100-G 

GENERAL RADIO COMPANY 

WEST CONCORD, MASSACHUSETTS, USA, 


4* Ca rrier Amplifier (Bridge Balance) 

The type 3 C 66 Carrier Amplifier plug-in unit is designed 
to be used with Tektronix 560-series osscilloscopes that will 
accept 3-series plug ip units. The unit is self contained and 
requires no external equipment other than the transducer. 
Excitation Voltage for the transducer is provided by the above 
unit. Suppressed -carrier amplitude modulation produced by unbalan- 
cing an A-C bridge with strain gages. A phase sensitive demodu- 
lator produces the proper direction of trace deflection. 


TECKTRONIX INC. BEAVERTON, OREGON (ESA). 


5- Strain Gages CA-10 (Bonded Type) 


Base Length 
Base Width 
Resistance 
Gage Factor 


- 10 mm 

- 6 mm 

- 111 . 0 + 6.5 ohms 

- 2.05 


Manufacturers Mahavir Optical & Scientific Works, 

44 » Civil lines, Roorkee, U.P. (India). 
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5.10 BESTJLT3 AND DISCUSSION 


The experimental results have been shown in Figs. (2.1 to 


2*10), The deflection at X = L/2, along at other points of "beam have 


been measured. At X = 1/2, L, l/5 results have been plotted against 
in Fig. (2.1, 2.1a, 2.1b). The measured values are in good agreement 
with theoretical values. Measured values are in between the values 
predicted by Modal Analysis and Energy Method. But the experimental 
values are in better agreement with Modal Analysis than Energy Method. 
In Fig. (2.2), Deflection has been plotted against time. The 
experimentally obtained deflection-time curve is also nearly 
sinuosidal in character. These traces have also been shown in plate g 
as obtained on the storage oscilloscope. At other points of the beam 
also similar observation have been obtained, though not shown. The 
deflection has been also plotted as magnification factor by dividing 
these values by Px/2K. 




The experimental value of strain at X = L /2, and X = L 
against v q has been plotted in Fig. (2.3, 2,5). These values also 
linearly increase with v Q . But experimental values aro higher than 
the analytically predicted values, at X = L/2, X = L/5. At X = L/2, 
the experimental values ore about 1 .4 times the value predicted by 
Modal Analysis and two times the value predicted by Energy Method. 
However at X = L, Fig. (2.5), the measured value of strain was O.7O 
times the value predicted by analytical method and two times the 
value predieted by Energy Method. At X = l/5 also the measured 
values are higher than the theoretical value. This difference 
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may "be due to hole in the experimental "beam, made for connection 
with the spring which are points of stress discontinuity* The 
strain-time curve has "been plotted in Fig. (2.4) and Fig. (2.6) 
for X = I1/2 and for X = 1 respectively and they are also shorn in 
plate 10. 

The impact force has been shown in Fig. (2.7,2.8). The 
maximum experimental value of impact force is less than the maximum 
value predicted by Modal Analysis and more than that predicted by 
Energy Method. The Acceleration and velocity time trace has been 
shown in plate (s and 7 ). 

For measuring the displacement and acceleration, calibration 
of the vertical scale of oscilloscope has been done and shown in 
Figs. 5*5 5*6) respectively. 

In case of strains at x - L/2, the difference between 
theoretical valaes and experimental valies may be due^ to three- 
dimensional effects of the central impact . 



TABLE 5.f 


LENGTH OP AIL BEAMS = 59.98 inches 
SPRING CONSTANT = 99.8 lh/in 2 


Beam Size 


No. 


D 

inch 


H 

inch 


lb/i 


xn 


E 

S i 5 

i 


mg 
^ ( 
lb. 


Hammer 
Q Radius 
0 inch 


BEAM I 

1.98 

0.495 

0.292 

50x1 o 6 

5. 75 

2 

BE M II 

1.98 

O.25O 

0.292 

30x1 o 6 

3.75 

2 


BEAM I 


1 

No. of ( 

| Natural Frequence = | 

Simply Supported 
Beam -without 
foundation 

Mode n ( 

| Free-Free 

| Simply Supported jj 

1 

Beam 

536 

on Elastic Foundation 

403 

173 

2 

1154 

785 

695 

3 

1827 

I6O7 

1550 

4 

3538 

2806 

2480 

5 

, 5270 

4364 

4320 

6 

7352 

6272 

6200 

7 

9783 

8531 

8448 

8 

12563 

11 138 

moo 

9 

15690 

14094 

14020 

10 

I92OO 

17399 

17300 


Coefficient of restitution for free-free heami ©n 
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CHAPTER 6 


CONCLUSIONS MD HECCMHENDATIOHS 

6.1 GENERAL CONCLUSIONS 

The response of "beam on elastic foundation has been 
analyzed by Modal Analysis and Energy Method. The entire accuracy 
of method using Modal Analysis depends upon the assumption of initial 
velocity distribution. The distribution assumed in this present analy- 
sis has not been truly represented in the experimental model, 
owing to the small length of the model beam. However under these 
limitation^ they are in good agreement. Impact load on elastic 
foundation is essentially a problem of dynamic plasticity , and 
permanent deformations are of considerable importance . Bounds can 
be. computed for such deformation using Energy Principles . 

Displacement Bounds within the elastic limit can be 
computed using Energy theorem and also by suitable inequalities 
and variational technique. The study of wave propagation along 
the beam on elastic foundation reveals the validity of assumptions 
usually made in Modal Analysis or Energy Method. 

The following conclusions have been made, based on the 
analytical and experimental results, obtained in the present 
investigation . 

(l) The Modal Analysis and Energy Method can be used 
to predict the response of beams on elastic foundation, subjected to 
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Impact loads. According to present study, Modal Analysis is 
preferable to Energy Method. 

( 2 ) For free-free "beam on elastic foundation, the rigid 
body translational mode end rotational modes modify the beam 
response. In the present study the deflection at the centre was more 
than at the free-free ends, however it is more a function of length 
of beam and it is also dependant rigid body modes and number 

of modes excited by. the impact load. 

( 3 ) The effect of foundation modulus is to increase the 
natural frequency of thebeam on elastic foundation. However this 
effect diminishes as number of mo des increases. 

( 4 ) For free-free beam, the first mode is predominating, 
and the deflection as well as strain decreases rapidly as the 
number of mode shape increases. However for simply supported beam, 
the first mode gives the major part of deflection, but for strains 
higher modes contribute more than the first mode. 

( 5 ) The external damping greatly reduces the deflection 
as well as strain. For £ < 1,0 the motion is oscillatory with 
diminishing amplitude. And for E > 1.0, the motion is no longer 
oscillatory being exponentially decreasing wi th. time. 

(6) The experimental values of deflection, impact force 
are in fairly good agreement with theoretical values predicted by 
Modal Analysis, and ' Energy Method. However in this study agreement 



is more pronounced with Modal Analysis than Energy Method. In the 
case of strains, the agreement is not that pronounced. 

(7) For a beam on elastic foundation subjected to dynamic 
load and impact load, where in initial velocity distribution is 
known, bounds on displacement can be computed using Energy Method* 
Within the elastic range, for a half sine curve the computed bound 
is only 12 fo higher than the exact value. These bounds are less for 
beams on elastic foundation than for beam without foundation. 
Assuming an elastic plastic moment curvature relationship, elastic 
contribution towards the bound is very small. So it can be safely 
taken as rigid-plastic, as far as the computation of these bounds 
outside the elastic range are concerned. The contribution from 
elastic range is more in beams without foundation than beams on 
elastic foundation. 

(8) Impact load caises permanent deformations on the 
beams. Without necessarily knowing the deformation process, the 
ultimate permanent deformation can be obtained by these bounds, 
computed outside the elastic range. With increase in foundation 
modulus these bounds decrease. 

(9) For a central impact by a mass m^ with a velocity v q , 
for a simply supported beam on elastic foundation, the computed 
bound deflection is only 1,5 times more than that by Modal Analysis. 
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(10) For a general class of loading q(x, t), displacement 
“bounds can be computed for a beam on elastic foundation as given 
in Eq* (3*5.28) andEa. (3*5*30) by using variational principle* 
For impact loads as given by Hertz, these bounds are about twice 
the value predicted by Modal Analysis and thrice the value 
predicted by Energy Method* 

(11) The beam on an elastic foundation is a dispersive 
medium for wave propagation. The flexural nave moves faster 

in a beam on elastic foundation than in a beam without foundation, 
where as energy propagates at a slower rate along a beam on 
elastic foundation than along a beam without foundation. 

6.2 SCOPE FOE FOE THEE BE SEARCH 

The following topics are listed for further investigations 

1) Impact loads on long beams on Elastic foundations. 

2) Experimental studies taking care* of Damping in the 
Foundation medium 

3 ) Experimental Methods for separating rigid body motion 
in the case of beams on elastic Foundations with free ends. 

4 ) Effect of shear Deformation androtatory inertia on 
the responses of beams on elastic foundations subjected to 
Impact loads. 
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